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Abstract 



We develop a method to solve the Bogoliubov de Gennes equation for 
superconductors self-consistently, using the recursion method. The method 
allows the pairing interaction to be either local or non-local corresponding to 
s and ci-wave superconductivity, respectively. Using this method we exam- 
ine the properties of various S — N and S — S interfaces. In particular we 
calculate the spatially varying density of states and order parameter for the 
following geometries (i) s-wave superconductor to normal metal, (ii) d-wave 
superconductor to normal metal, (iii) ci-wave superconductor to s-wave super- 
conductor. We show that the density of states at the interface has a complex 
structure including the effects of normal surface Friedel oscillations, the spa- 
tially varying gap and Andeev states within the gap, and the subtle effects 
associated with the interplay of the gap and the normal van Hove peaks in 
the density of states. In the case of bulk <i-wave superconductors the surface 
leads to mixing of different order parameter symmetries near the interface 

and substantial local filling in of the gap. 
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I. INTRODUCTION. 



Interfaces in superconductors, especially high temperature superconductors, are of con- 
siderable interest for both fundamental physics and for applications. Many experiments 
have used both single electron tunnelling and Josephson effects as probes of the energy gap 
and order parameter symmetry in the cuprates. For example Wollman et al. [T|] and Sun 
et al. constructed SQUID devices consisting of junctions between YBa2Cus07 (YBCO) 
and Pb, while Tsuei et al. || constructed superconducting rings consisting of YBCO thin 
films with two or three grain boundary junctions. Theoretical analysis of experiments such 
as these relies to a large extent on macroscopic symmetry arguments and not on the micro- 
scopic details of the actual interfaces. However in some cases the microscopic physics at the 
interface can be an important factor in understanding the experimental results. For exam- 
ple, if mixing of different order parameter symmetries occurs at the interface (because the 
interface breaks the bulk tetragonal or orthorhombic symmetry) the extent of such mixing 
can only be determined from microscopic calculations. Similarly, suppression of the order 
parameter (either <i-wave or s-wave) near an S — N interface can lead to significant local 
density of states within the bulk energy gap, and this can complicate the analysis of single 
electron tunnelling spectra ||. The microscopic physics of surfaces and interfaces of high 
T c superconductors are especially interesting because of the short coherence length, and the 
probable <i-wave gap function. 

In the past few years there have been a number of microscopic calculations of surfaces 
and interfaces in superconducting systems with a <i-wave order parameter. Most of the 
theoretical results have been obtained using tunnelling theory, or Andreev's approximation 
fl5HTTl in which the tunnelling barrier and the order parameter are not found self-consistently. 
For tunnel junctions these approximations may be adequate, but we show below that self- 
consistency has significant effects for interfaces with direct contact between the constituents. 
Self-consistent properties of interfaces have previously been computed using the Eilenberger 
equations [f[^-|PT||, which are an approximation to the Bogoliubov de Gennes equation. These 
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self-consistent solutions to the Eilenberger equations have shown some interesting effects 
which have only arisen when the order parameter is calculated in a self-consistent manner. 

In this paper we aim to show how to calculate self-consistent properties of superconduct- 
ing interfaces by directly solving the Bogoliubov de Gennes equations. This approach has 
the advantage that self-consistency can be fully incorporated and also there are no need for 
the further approximations of the Eilenberger method. Our approach makes use of the recur- 



sion method [|18j to solve the Bogoliubov de Gennes equation on an arbitrary tight binding 
lattice. Previously the recursion method has been used to examine the effects of disorder 
in s-wave []l9| , pQ| and ci-wave |2l]] superconductors, and to determine the core structure of 



vortices and explain the origin of de Haas van Alphen oscillations in the superconducting 
state ||22|| . In particular Litak, Miller and Gyorffy have given a detailed description 
of the application of the recursion method to local interactions, corresponding to s-wave 
superconductors. Here we extend this to the case of non-local interactions necessary to ob- 
tain d-wscve superconductivity and apply the method to various interfaces of s and <i-wave 
superconductors. 

The method of performing our self-consistent calculations will be described in section [IT]. 
Here we introduce the Bogoliubov de Gennes equation with a general interaction, Uy, and 
demonstrate how this general Hamiltonian can be solved self-consistently using the recursion 



method ||T8| , |20| . A necessary and non-trivial step in the calculation, as described below, 
involves finding the density of states accurately by extrapolation of the recursion method 
continued fraction. 



In Sec. m we proceed to apply the method to several different problems. Firstly various 
test calculations are described, including: the local density of states for a uniform system 
with no interactions, the local density of states for a system with a local attractive interaction 
(local s-wave superconducting order parameter) and a system with a non-local attractive 
interaction. We show that for this non-local interaction there are two possible solutions, these 
solutions being a non-local s-wave superconducting order parameter (extended s-wave) and 
a non-local c?-wave superconducting order parameter. 



Having tested the method on uniform systems we present our self-consistent solutions 
for the interface between two different materials. We will consider three different interfaces. 
First we consider a normal metal to s-wave superconductor, N — S s , interface where the 
pairing interaction is zero in the normal region (N) and purely local in the superconducting 
region (S s ). Then we will consider a rf-wave to normal metal, S d — N, interface where again 
the interaction in the normal region is zero and there is a non-local attractive interaction in 
the superconducting region (S d ). Finally a study of a S d — S s interface will be described. 
These three calculations enable us to make a comparative study of how the local density of 
states and order parameter changes as a function of position across the different types of 
interface. 



II. THEORY/MODEL. 



A. The Bogoliubov de Gennes Equation. 



The Bogoliubov de Gennes equation on a tight binding square lattice has the form 
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where w™ and vf are the particle and hole amplitudes, on site i, associated with an eigenen- 
ergy E n and where is the (possibly non-local) pairing potential or gap function. 

In the fully self-consistent Bogoliubov de Gennes equation the normal state Hamiltonian 
is given by 



= (Uj + ^Uijnij){l - 5ij) + (cj - n + \Vini)bij 



(3) 



where \x is the chemical potential, e« is the normal on site energy of site % and tij is the 
hopping integral between site i and site j, for the rest of this paper is non zero for 



nearest neighbours only. The on-site and off-site interaction terms \Uunu and ~£/y7iy are 
the Hartree-Fock potentials corresponding to the on-site interaction Ui and the non-local 
interaction [/y. The charge density entering the Hartree-Fock terms riy is given by 

mj = £<*U*> = 2E(W*«"/(^) + «?(«")*(! - f( E n))- (4) 

a n 

Similarly the pairing potentials are defined as 

Ay = -UijFij (5) 

where the anomalous density is 

Fa = = £««)*(i - /(^)) - K n )*K)/(K). (6) 

n 

In equations (|4]) and (0) the sums only consider terms E n up to the condensate chemical 
potential (/i). 

A solution to the above system of equations will be fully self-consistent provided that 
both the normal (fTy/iy) and anomalous (Ay) potentials are determined consistently with 
the corresponding densities ny and Fy via Eqs. ^ and [5| Note that the normal Hartree- 
Fock terms £/yny play an important role and cannot be neglected. For on-site interactions 
these terms correspond to position dependent shifts in the on-site energy, while for non-local 
interactions these terms renormalise the hopping £y leading to position dependent changes 
in the electronic bandwidth. 

Figure 1 illustrates the geometry corresponding to this system of equations. The tight- 
binding lattice has nearest neighbour hopping interactions {tij), as well as a coupling between 
particle and hole space, via a superconducting order parameter (Ay). If the interactions 
are purely on-site (Ui) attractions then the pairing potential will be purely local (A&), 
corresponding to the dashed line in Fig. 1. On the other hand when the interaction is 
non-local (Z7y, % ^ j) the pairing potential Ay will also be non-local, as illustrated by the 
solid lines in Fig. 1. For computational convenience we limit both the hopping and non-local 
interaction to nearest neighbours distances. We also need to specify over what energy range 
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the interaction has an effect, and as in BCS p3 |, we will assume that it only acts over a 
small energy range centred on the Fermi energy, /z ± E c . 



B. The Recursion Method. 



The method we have adopted to solve the above system of equations is the recursion 



method [181. This method allows us to calculate the electronic Green's functions 



G aa ,(i,j,E) = (ia\ 



1 
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where the indices i and j denote sites, while a and a' represent the particle or hole degree 
of freedom on each site. We denote particle degrees of freedom by a = + and hole degrees 
of freedom by a = — . For example j, E) represents the Greens function between the 

particle degree of freedom on site i and the hole degree of freedom on site j. 

To compute the Green's functions (|7|) we can closely follow the method described by 
Litak, Miller and Gyorffy [^UJ for the special case of a local interaction (Uij = UuSij). Using 
their method we can transform the Hamiltonian to a block tridiagonal form 
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where a n and b n are 2x2 matrices. Given this form for (ia\El — H\ja') and expressing 
the Green's function as 



G ao ,(i,j,E) = (ia\(El-H)- l \ja') 



(9) 



the Greens functions above can be evaluated as a matrix continued fraction so that 
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G{i, j, E) = (El -a -b{ (El -a t -b\ (El - a 2 - bj (El - a 3 - . . .) 1 b 3 ) ' 6 2 ) * b^j 

(10) 

where 



G(i,j,E) 



* G aa (i,i,E) G aa ,(i,j,E)\ 



(11) 



yG a/a (j,i,E) G a < a >(j,j,E) j 
Within equations (§) and (|I0|) we have a formally exact representation of the Green's 
functions. However in general both the tridiagonal representation of the Hamiltonian, and 
the matrix continued fraction @ will be infinite. In practice one can only calculate a finite 
number of terms in the continued fraction exactly. In the terminology of the recursion 



method it is necessary to terminate the continued fraction p8|j20| , |2^ -t28| . 

If we were to calculate up to and including a n and b n and then simply set subsequent 
coefficients to zero then the Green's function would have In poles along the real axis. The 
density of states would then correspond to a set of 2n delta functions. Integrated quantities 
such as the the densities n^ and Fy could depend strongly on n, especially since only a few 
of the 2n delta functions would be within the relevant energy range within the BCS cut off, 
E c . In order to obtain accurate results it would be necessary to calculate a large number of 
exact levels, which would be expensive in terms of both computer time and memory. 

As a more efficient alternative we choose to terminate the continued fraction using the 



extrapolation method, as used previously by Litak, Miller and Gyorffy f20|| . We calculate 



the values for a n and b n exactly up to the first m coefficients using the recursion method. 
Then, noting the fact that the elements of the matrices a n and b n vary in a predictable 
manner ||20|| , we extrapolate the elements of the matrices for a further k iterations, where k 



is usually very much greater than m. This enables us to compute the various densities of 
states, and the charge densities n^ and Fij accurately with relatively little computer time 
and memory. 

In terms of the Green's functions G aa /(i, j, E) the pairing and normal Hartree-Fock 
potentials A^- and hUijUij are given by 
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**i = TT U * I ° (G + -(h3,E + ir])-G + „(i,j,E -i V ))(l- f(E))dE (12) 

ZTT J-Ec 

and 

\Uijnij = ^Uij £J (G ++ (i,j, E + irj)- G + + (iJ, E - i V ))f(E)dE. (13) 

where rj is a small positive number. 

To obtain the above equations we have used the fact that 
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are the eigenvectors of equation ([[]) with eigenvalues i? n and — E n has been used. Also note 
that the integrals in equations (0) and (0) are bounded by the cut-off E c , corresponding 
to the energy dependent interaction 

I -\U\ for \E-u\< E c 
U lJ {E) = \ i (15) 

! for \E-y\> E c 

as in BCS theory. Our cut-off E c can correspond to the BCS cut-off Hud arising from 
retardation of the electron-phonon interaction, or any other energy scale cut-off for the 
interaction which may be applicable for high temperature superconductors 



C. Achieving Self-consistency. 

Using the above methods to calculate Ay and UyTiij we need to achieve a fully self- 
consistent solution. Firstly we make use of any symmetries in the system in order to minimise 
the number of calculations which are necessary. For example on an infinite square lattice 
with no variation in any of the potentials only one independent site needs to be calculated 
since this site can be mapped onto all of the other sites. Secondly, once we have decided 
which sites need to be calculated self-consistently A^ and can be calculated for those 
sites, remembering that on a square lattice each site will have four nearest neighbours. This 
implies that in general we will have to calculate equation nine different Green's functions in 
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order to calculate A^ and n^. This can be seen by considering site i in figure 1 and noting 
that we need to calculate the Greens functions shown in table I, depending on whether the 
interaction is purely local, purely non-local, or both local and non-local. Having calculated 
the appropriate Green's functions new values for A y - and n y - can be calculated, which we will 
denote as A^ and n$ . Inserting these into the Hamiltonian and repeating the calculation 



(2) (2) 

of then Green's functions leads to a new set AJ - and n\ ■ and so on. We repeat this iteration 



for all i and j until 



and 
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< 0.001 



< 0.001. 



(16) 



(17) 



Since A y - and n y - can be complex we need to also check for convergence in their asso- 
ciated phases. We do this and find that convergence in the phase gradient of the complex 
parameters is much more rapid than the convergence in the magnitude. 



III. NUMERICAL RESULTS. 
A. Uniform Systems. 

As a first test of the above methods let us examine a bulk superconductor, corresponding 
to an infinite 2-d square lattice with either local or non-local attraction. These examples 
will show how well such quantities as the local particle density of states can be calculated 
and how the extrapolation of the elements of the matrices a n and b n is performed. 

The quantity of interest is the local particle density of states, which can be calculated 
from the following expression 

N t (E) = ^(G + + (i, i,E + i V ) - G + + (i, i,E — i V )). (18) 

Consider first a non-interacting system where = 0, \i = 0, = 0, and = 1 for 
nearest neighbours and zero everywhere else. Figure 2 shows the local particle density of 
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states for a calculation where the number of exact continued fraction levels was m = 50 and 
the elements of a n and b n were extrapolated for 2000 more values. For this calculation the 
convergence parameter r\ was chosen as i] = 0.02. Fig. 2 shows that using the extrapolation 
method the central logarithmic van Hove singularity and the sharp band edges can be 
resolved very well. Figure 3 shows the first 100 continued fraction coefficients where the 
first 50 are calculated directly using the recursion method and the rest are the extrapolated 
values. It is clear from the figure that the oscillations in Wj 1 ^ still persist after the first 50 
continued fraction levels. In fact these oscillations die off slowly as 1/n and it is critical 
to include them correctly. From figure 3 it is clear that the first 50 levels provide enough 
information about the decaying oscillation that the Wy^ can be extrapolated quite easily. 

Having considered a system where the interaction is zero the next step is to consider 
systems where the interaction is uniform and finite. For such systems the local particle 
density of states can be calculated, for different types of interaction. In figure 4 we have 
plotted two different local particle densities of states for a local interaction = — 2.55^ 
(dashed line) and Uy = — 2.5(1 — 5^) for nearest neighbours (solid line). In each case E c = 4 
and tij = 1 (for nearest neighbours) and all other parameters were set to zero throughout 
the lattice. 

The dashed line in figure 4 clearly shows the energy gap at the Fermi energy, characteristic 
of s-wave superconductivity. The van Hove peak in the density of states is also very clearly 
resolved. The solid line in Fig. 4 shows the local particle density of states going to zero at the 
Fermi energy, in a manner which is typical of the local particle density of states for a rf-wave 
superconductor. In this case of the non-local interaction the order parameter changes sign 
we as rotate by n/2 around a site, ie. in reference to figure 1 = — Ajj 2 , A^ 2 = — A^ and 
Aij A = — Ajj 3 . The way we have performed the calculation is to keep the Fermi energies the 
same in the two calculations but change, in the case of the local interaction (dashed line), 
the density and, in the case of the non-local interaction (solid line), the width of the band 
self-consistently. This has the effect of moving the system away from half filling in the case 
of the local interaction, and, broadening the band in the case of the non-local interaction, 
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because of the local and non-local Hartree-Fock terms in the Hamiltonian Uanu and Uijfiij 
respectively. 

At this point one should note that in the case of a non-local interaction one can as well 
as having a <i-wave self-consistent solution to the Bogoliubov de Gennes equation, it is also 
possible to obtain an extended s-wave solution, ie. = Ajj 2 , A^ 2 = A^ and Ay 4 = Aj 3 - 3 . 
However we find that such solutions are less stable than the rf-wave solutions, this is only 
true at or near half filling of the band. 

To obtain the results shown in figure 4 we have again calculated 50 levels of the recursion 
method exactly and then extrapolated for further 2000 levels. This can easily be done 
because the elements of a n and b n vary in a predictable manner, as has already been seen 
for the case without interactions. 

B. Interfaces. 

Having considered systems where the interactions remain uniform throughout the struc- 
ture the next step is to consider systems which contain interaction strengths which vary in 
real space. The most simple case one can conceive for this scenario is an interface. We will 
simply model the interface by allowing the interaction to change in a step like manner. 

We will consider three separate situations, N — S s , S d — N and S d — S s . In the normal 
region we shall set Uij = 0, hence the order parameter in this region will be zero, (but one 
should note that this does not imply that is zero). Before we look at the numerical results 
it is worthwhile considering what one may expect to find. In the case of a local interaction 
the results are well documented, i.e. the magnitude of the superconducting order parameter 
reaches a maximum at the bulk value a few coherence lengths in the superconducting region 
away from the normal interface. In the case of the non-local interaction we would also 
expect the amplitude of the superconducting order parameter to reach a maximum several 
coherence lengths away from the interface, but the problem of how to define the magnitude 
of the superconducting order parameter now arises. Going back to figure 1 we can see that 
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for each site % there are five Ay's, so hence for each site we can define five order parameters 
per site. We can also combine these different order parameters on each site in the following 
manner 

| A ( S (/ocaO)| = | A .| ( 1Q ) 

|Af I = \\^in ~ + Ay, - Ay 4 | (20) 

^(non-local^ = ^ + ^ + ^ + ^ (21) 

so that each equation defines a different type of symmetry for that site. Since the systems 
we are interested in change in the x-direction only it is possible, when one is considering the 
properties of that interface, to look along one line of sites in the x-direction and note that 
for any other y coordinate the properties of the system are the same, so Aj — > A(x). 

Having defined all the quantities of interest the next step is to specify some of the systems 
of interest. The three systems we are going to consider are as already pointed out N — S s , 
S d — N and S d — S s , to set up these systems we used the parameters shown in table II. 

Figures 5(a-c) plot the three main symmetry components of the order parameter 
|AW Ioc °0)(a;)| (dashed line), \A^(x)\ (circles) and \A (s(non ~ local »(x)\ (solid line) for the three 
different geometries N - S s (figure 5(a)), S d -N (figure 5(b)) and S d -S s (figure 5(c)). The 
interface corresponds to x — 100 on the figures. Figure 5(a) shows, as expected, that the 
s-wave order parameter, \A^ local ^(x)\, simply rises over a coherence length to a maximum 
at the bulk superconducting order parameter. Because the interaction is purely on-site in 
Fig. 5(a) \A^(x)\ = \A( s ( non - local » (x)\ = 0. 

In figure 5(b) we see that for the <i-wave to normal metal interface |A^(x)| also drops 
to zero at the interface. However, unlike the s-wave case, it does not simply drop to zero 
smoothly but has a sharp peak structure right at the interface. The origin of this peak is 
explained by looking at the extended s-wave component, | /\( s ( non - l ° cal )) ( x )\ (solid line in Fig. 
5(b)). We see that the extended s-wave gap function is finite near the interface. This is due to 
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he fact that the order parameter varies near the interface and hence (x) ^ Aj 3 (x), making 
the values of | A( s ( no ™ _k,ca '))(:c)|. This is emphasised in figure 5(d) where A J - 3 (x) — A jl (x) 
is plotted, from this graph one can see that the peak in |A^(x)|, in figure 5(b), near the 
interface is due to the component in the the x direction. 

Figure 5(c) shows the ci-wave to s-wave superconductor interface. Again we can see that 
the extended s-wave component 

\ A (s(non-local)) 

(x) | is non-zero at the interface, even though 
it is zero in the bulk on both sides, and that this leads to sharp features in both the local 
s-wave and <i-wave order parameters near the interface. 

Having seen how the profiles of the superconducting order parameters are affected by 
the proximity of different materials, we now look at how the local particle density of states 
changes as we move across the various interfaces. Figures 6, 7 and 8 are contour plots of the 
local particle densities of states for the three interfaces of interest. Figure 6 shows a contour 
plot for the N — S s interface. Looking at this plot one can see that as we move across the 
interface, at x — 100, the superconducting gap opens up within a couple of atomic sites. 
On the normal metal side, for x < 100, the van Hove singularity in the centre of the band 
can be clearly seen, but as we move into the superconducting region the band edges are 
shifted (due to the Hartree-Fock potential term) and the superconducting gap opens up at 
E — 0. In the superconducting region the van-Hove singularity is shifted away from E — 
as can also be seen in figure 3 (dashed line). Due to the mismatch in the band edges we 
see oscillations in the local particle density of states near the band edges; these are simply 
Friedel oscillations. 

Figure 7 shows a similar contour plot of the local particle density of states for an S d — N 
interface. Again we can clearly see the gap in the superconducting region and the van Hove 
singularity in the normal region. In this system the Hartree-Fock potential term leads to an 
increase in overall band width on the d-wave side. Again since the band edges do not match 
up we see Friedel oscillations in the local particle density of states near the band edges. 

Finally in figure 8 we have plotted the local particle density of states as we move across 
the S d — S s interface. This plot has many interesting features, the first to note is that again 
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due to the mismatch in the band edges oscillations appear in the local particle density of 
states. Secondly for x < 100 (Sd region) the density of states gradually goes to zero at E — 
(typical of d-wave superconductivity (see figure 4 (solid line))). Whereas for the S s region 
the local particle density of states drops to zero very sharply. The main points of interest 
is what happens at the interface itself. In plane of the interface there are states in the gap, 
as both the d-wave and s-wave order parameters are suppressed. At x = 100 there are two 
peaks in the density of states just above and below E = 0, which as we move further into 
the S s region are shifted to become the BCS density of states singularities just above and 
below the superconducting gap. Note that the parameters for the calculation in Fig. 8 were 
chosen so that |A^| >> /\( s ( local ))\ as would be the case for a YBCO-Pb junction such as 
those used by Wollman et al. |J. 



IV. CONCLUSIONS. 

In this paper we have shown how it is possible to perform self-consistent calculations 
of the Bogoliubov de Gennes equation, using the recursion method. This method has the 
advantage of being an order N method and hence allows us to tackle problems with a 
relatively small amount of computational effort. A key to obtaining accurate densities of 
states with relatively little computational effort is the extrapolation procedure we have used 
to terminate the matrix continued fraction. Our method is fully self-consistent, including 
both self-consistency in the order parameter and in the normal Hartree-Fock potentials. 
As we have shown these normal potentials make significant contributions by shifting or 
widening the density of states in a spatially dependent manner. Our method can deal with 
both local attractive interactions, corresponding to local s-wave superconductivity, or non- 
local interactions corresponding to <i-wave or extended s-wave pairing. In our system we 
found that the (i-wave is more stable. 

As a first application of the method we examined three simple interfaces, corresponding 
to an s-wave S — N junction, a d-wave S — N junction and an s-wave to <i-wave S — S 
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junction. The numerical results show a number of interesting features, including a non- 
monotonic variation of the order parameters near the interface, a surface layer of extended 
s-wave pairing (even though it is not stable in the bulk), and subtle effects of the self- 
consistent Hatree-Fock terms in the Bogoliubov de Gennes Hamiltonian leading to Friedel 
oscillations and spatially dependent shifts in the van Hove singularities near the interfaces, 
as highlighted by the contour plot in figure 8. 

In future we hope to apply our method to more complex interfacial phenomena in super- 
conductors, such as junctions carrying supercurrent (e.g to look for tt-j unctions), supercon- 
ducting twin boundaries and grain boundary junctions. Our methods can also be applied 
to many other problems in superconductivity such as the structure of vortex cores in s or 
ci-wave superconductors, the effects of impurities and so on. 
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FIGURES 

FIG. 1. This is a schematic diagram of a tight binding lattice, with particle and hole degrees 
of freedom, Ajj couples particles on site i to holes on site j. The difference between local and 
non-local pairing is highlighted by the dashed (local pairing) and solid (non-local pairing) lines. 

FIG. 2. The local particle density of states for a 2-D tight binding lattice with no interactions 
(Uij = 0). For this system tij = 1 for nearest neighbours only, = and = 0. 

FIG. 3. A plot of the real part of b(n) 11 for the same system which was used to calculate the 
local particle density of states in figure 2. 

FIG. 4. Two plots of the local particle density of states for, a local interaction (Uijdij = —2.5, 
Uij(l — 5ij) = (dashed line)) and a non-local interaction (Uijdij = 0, Uij(l — 5{j) = —2.5 (solid 
line)), all the other parameters are equal to those used to obtain figure 2. 

FIG. 5. Figures 5(a-c) plot the profiles of different symmetries of the superconducting order 
parameter (|A( s ( non - /oca/ ))( x )| (solid line), |A( s (' oca/ ))( x )| (dashed line) and |A^(x)| (circles)) for 
different interfaces. Figures 5(a), (b) and (c) are for N — S s , S d — N and S d — S s interfaces 
respectively. Figure 5(d) plots Aj 1 (x) — Aj 3 (x) for the N — S d interface. The parameters used to 
obtain the figures are given in table II. 

FIG. 6. This is a contour plot of the local particle density of states as a function of position, as 
one move across the N — S s interface. The steps in the contour plot are in units of 0.4, ie. white 
represents N(E) < 0.4 and black represents N(E) > 0.32. The parameters used to obtain this 
graph are given in table II. 

FIG. 7. This is a contour plot of the local particle density of states as a function of position, as 
one move across the S d — N interface. The steps in the contour plot are in units of 0.4, ie. white 
represents N(E) < 0.4 and black represents N(E) > 0.32. The parameters used to obtain this 
graph are given in table II. 
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FIG. 8. This is a contour plot of the local particle density of states as a function of position, as 
one move across the S d — S s interface. The steps in the contour plot are in units of 0.4, ie. white 
represents N(E) < 0.4 and black represents N(E) > 0.32. The parameters used to obtain this 
graph are given in table II. 
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TABLES 



Interaction Type 


G+-(i,i,E) 


G+±(i,h,E) 


G+±(i,j 2 ,E) 


G+±(i,j 3 ,E) 


G+±(i,jt,E) 


U u 


Y 


N 


N 


N 


N 


Uij{\ - S^) 


N 


Y 


Y 


Y 


Y 


Ua + Uijil-Sij) 


Y 


Y 


Y 


Y 


Y 



TABLE I. This table shows which Greens functions need to be calculated for systems with 
interactions which are local, Uu, non-local, Uij(l — 5 — ij), or both. The site labels correspond to 
the notation of Fig. 1. 





x < 100 


x > 100 


System 


Uu 


Uij(l - Sij) 


E c 


Uu 


Uij(l-5ij) 


E c 


N-S s 











-2.5 





4.0 


S d -N 





-3.5 


4.0 











S d_ s s 





-3.5 


4.0 


-2.5 





4.0 



TABLE II. This table defines how the interactions vary in real space for the three different 
interfaces. All energies are given in units where the nearest neighbour hopping tij = 1. Also fx = 
everywhere and T = 0.01. 
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